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c/3 ' We consider a system of two species of bosons of equal mass, with in- 

teractions C/ a (|x|) and C/ X (|x|) for bosons of the same and different species 

respectively. We present a rigorous proof-valid when the Hamiltonian does 

not include a species switching term-showing that, when C/^dxl) > U a (\x 

the ground state is fully "polarized" (consists of atoms of one kind only). In 



the unpolarized phase the low energy excitation spectrum corresponds to two 
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O ■ polarization instability is signaled by the vanishing of the velocity of the odd 



linearly dispersing modes that are even and odd under species exchange. The 



modes. 
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The experimental observation of Bose-Einstein condensation (BEC) in dilute atomic 
systems [|IJ has triggered a very intense theoretical activity ||. Attention has broadened 
to include condensates with internal degrees of freedom, or multi-species BEC, which were 
realized for trapped rubidium || and sodium [[|. Early theoretical analysis of m-species 
condensates focuses on the large m limit for Hamiltonians invariant under U(m) transfor- 
mations II, superfluid helim mixtures || and spin-polarized hydrogen 0. In the context of 
BEC much of the theoretical attention concentrates on spinor condensates || in which the 
internal degrees of freedom correspond to the different Zeeman states of the of a particular 
hyperfme manifold such as the / = 1 manifold in sodium. In these cases the two-body inter- 
action is invariant under under rotation in species space. On the other hand, for rubidium 
one has two non-degenerate internal states corresponding to two different manyfolds, and 
the two-body interaction is not invariant under SU(2) rotations in species space 0. The 
ground state and excitation spectrum of a two-species (a and b) condensate of this kind 
was studied in mean field in Ref. |10|| , where it was shown that the quasiparticle energy 
can become imaginary, signaling an instability when the inter-species repulsion U x is larger 
than the intra-species repulsion U a . This kind of treatment follows-as do the majority of 
theoretical approaches to BEC-the Gross-Pitaevskii (GP) mean field equation [IT]. On the 



other hand, rigorous results for BEC (and many body problems in general) are scarce and at 
the same time useful in providing control for approximate solutions. With this motivation, 
in this paper we consider a two-species system of interacting bosons, and show rigorously 
that the instability mentioned above corresponds to the tendency of the system to "polar- 
ize", the true ground state consisting of only one species of bosons when U x > U a . We also 
discuss the low-energy excitation spectrum and show that the instability is signaled by a 
divergence of the "compressibility" associated with exchanging particles from one species to 
the other at fixed total particle number. 

The Hamiltonian for our two-species system (a and b) of bosons of equal mass m is given 
by (h = 1) 



H = K + U a + U b + U ab , (1) 

with 

K = E ^ K W a (k) + WW^(k)] (2) 

Ua + Uh = \S d * xdZ v ua (\* - yl) M x K(y) + P6(x)p 6 (y)] (3) 

^ = |rf 3 xrf 3 y^(|x-y|)p a (x)p 6 (y). (4) 

In the above equations X is the kinetic term, and the terms (U a + £/&) and U ab corre- 
spond to the interaction between bosons of the same and different species respectively. The 
operators ^(x) destroy a boson of species i (i = a,b) at position x, and obey the follow- 
ing commutation relations [ , 0J(x),'0-(x')] = SijS(jc — x'). Also, Pi(x) = ■0j(x)^j(x), and 
■0i(k) = V -1 / 2 /d 3 x exp(ik • x)^ i (x) with V the total volume. 

Using a variational argument we will prove that, for potentials satisfying C/^dxl) > 
£/ a (|x|) the ground state of the above Hamiltonian is completely polarized. By "polarized" 
we mean that either of the two situations is realized: \N a ) = 0,(N b ) = N; or \N a ) — 
N, \Nb) — 0, with N the total particle number and iV, = / d 3 x Pi(x). 

We note that, since the masses are the same for both species, the kinetic term of Hamil- 
tonian (|J) commutes with the "rising" operator defined as: 

O r = /tfWi(x)V 6 (x) = E^(k)^(k) (5) 

J k 

which conserves the total number of particles but converts particles of type b into particles 

of type a. More specifically, for the case of different masses for each species we have that 

[^^] = fe-i)E k ^t(k)A(k). 

We have in mind the alkali atoms, which have a hard core interaction, meaning that the 
exact wave functions vanish when the coordinates of two atoms of either species coincide. 

Since the Hamiltonian conserves the particle number for each species, we can start with 
the normalized ground state wave function |\l/o(^a, iV&)) in the subspace of N a (N b ) particles 
of species a(b). The unpolarized situation corresponds to N a = N b = N/2. Now let us 



consider the normalized variational wave function \^ v (N a + N b ,0)) obtained by the action 
of the rising operator on \^ (N a , N b )) N b times: 

\* v (N a + N b ,0)) = -^=(O n ) Nb \MN a ,N b )). (6) 

Note that *$? v represents a completely polarized wave function, with particles of species 
a only. Since [On, K] — 0, we have 

(V v (N a + N b ,0)\K\V v (N a + N b ,0)} = (y (N a ,N b )\K\y (N a ,N b )}, (7) 

meaning that the completely polarized variational wave function and the ground state of the 
subspace (N a , N b ) have the same expectation value of the kinetic energy We stress that the 
function \ty v (N a + N b ,0)) as defined in Eq. @ is normalized only because the exact wave 
function \ty (N a , N b )) vanishes when any two coordinates coincide. Otherwise we would have 
to worry about permutation factors whenever coordinates coincide. 

In order to compute the change in the potential energy we write the expectation values 
of U a and U x as an integral over all the multiparticle configurations V = (r a , T b ) with 
coordinates {X}r = {~x.i t r a , ■ ■ ■ ,XN a ,r a ',xi,r b , ■ ■ ■ ,*-N b ,r b }- For each configuration T let us 
regard the particle coordinates {X} r as nodes of a graph. There are N a (N a — l)/2 and 
N b (N b — l)/2 edges connecting pairs of particles of species a and b respectively, and N a N b 
edges connecting a particles of different species. The contribution to the expectation value of 
the potential energy U = (U a ) + (U b ) + (U ab ) from this configuration is a function of the length 
of the edges of the graph, which can be classified in three sets: {^ir a > ' ' ' ^N a (N a -i)/2r a }y 

Ui,iV ' ' 'i&N b (N b -i)/2,r b h and W!n ■ ■ -^aNt^h where ^r„ is one of the possible lenghts 
|xfc,r a — x z,rJ, etc. The potential energy is therefore given by 



U = Jd 3N x r \y (x r ) 



'N a (N a -l)/2 N b (N b -l)/2 N a N b 

i=l i=l i=l 



(8) 



with \&o(Xr) the ground state wave-function in first quantization. 

In the variational wave function all the edges of type b and ab are both changed to edges 
of type a. Therefore the contribution to the potential energy of each edge of configuration T 



changes according to £7 a (^ r J -► £/ a (q, r J, U x (fifr) -► ^ a (C r )- If we cal1 ^ the expectation 
value of the potential energy in the variational wave function, and E v the variational wave 
function, we obtain that AE = E v — Eq = U v — Uq is given by 

. N a N b 

AE = / d N X v |vl>o(Xr)| 2 £ \U a {^r) ~ ^fr)] ■ (9) 

We see that when the interactions satisfy C/^dxl) > t/ a (|x|) for all values of x, AE 
is a sum of negative terms. In general the potentials have a repulsive short-range term 
and a long-range attractive tail. The condition for validity of our proof is that there are 
no "crossings" of the potentials [/ x (|x|) and C/ a (|x|) as a function of the coordinate. The 
simplest approximation will be to take the attractive components of both potentials as 
equivalent (of the form —Cq/R 6 ) and differing short range components. This dependence is 



consistent with calculations for ultra-cold Na collisions [12| . 

Since the completely polarized variational wave function has lower energy than the true 
ground state of the partially polarized subspace, the gound state will be completely polarized. 
It is also evident that the proof is also valid in the case different intra-species interactions 
U a (x) and U b (x). As long as U x > U a ,U b , the ground state is polarized with particles of 
type a {b) when £/ a (|x|) < £/ b (|x|) ((t/ 6 (|x|) < £/ a (|x|)). This means that the difference 
C/ 6 (|x|) — £/ a (|x|) plays the role of a "symmetry breaking field". Also, note that the above 
proof is also valid in the presence of an external potential U e (x) that is equal for both species, 
meaning that the polarization transition also occurs for trapped atoms. 

We note that the bosonic nature of the particles is crucial for establishing our rigorous 
proof. If the atoms a and b were fermions, the variational argument ceases to be valid: since 
the Hamiltonian conserves species, atoms belonging to different species can be considered 
distinguishable, and the wave function does not change sign if we exchange any two atoms 
a and b following a path P. If we convert an atom a to an atom b, the wave function 
has to change sign under the particle exchange following the same path P, implying that 
we have to introduce an additional node in the wave function. Formally this means that 
O R \V {N a ,N b )) = for fermions. 



We now discuss the low energy excitation spectrum, which in the symmetric case (U a = 
U b ) can be computed using longitudinal sum rules [[U||l4]]. In the unpolarized case the 
excitations correspond to two phonon branches with wave functions p k = (p k ± Pk)|0). The 
operators p k = V~ x l 2 J <i 3 o;e ik ' x p Xj ; are the Fourier transforms of the density operators for 
each species. In other words, since the Hamiltonian is invariant under exchange of species, 
the excitations are either even or odd in the species index. 

The excitation energies are given by 

k - w&m ' (10) 

where we have shifted the origin of energies (H —* H — Eq\ H\0) = 0). We will first show 
that (0|p± k #p£|0) = nk 2 /2m, with n = N/V. Consider 

(0\p\Hp b k \0) = i(0| {[p\,H]pi + p\[plH]} |0> (11) 

= |r E k ■ q(o| K(q' - k)^(q')4(q + k)A(q) 

-$(q + k)^(q)^(q' - k)^(q')l |0>, (12) 



with Eq. ([I4) following from the direct evaluation of the commutators in Eq. (|Tl|). 
Since the ground state is symmetric under species exchange, and since the operators b 
and a commute, Eq. flTJD implies that (0|p^ k #p k |0) = (0|p*l k iJp£|0) = 0. On the other 
hand, if in Eq. (|T2|) we replace b with a we obtain the usual sum rule (0|p" k ifp k |0) = 
(h k 2 /2m)V~ 1 (0\ 5Z q Va(q)Va(q)|0)) which implies that for the even and odd modes the 
exact /-sum rule applies |H|: (0| p_ k Hp k |0) = nk 2 /2m, or, which is equivalent: 

duJUjS ± (k,uj) = , (13) 

2m 

with 

$*(*, w) = n- 1 ^ I (HPk |0> Moo - uo u ) (14) 

1/ 

the dynamic structure factor of the even and odd modes. The excitation energies are there- 
fore given by 



^ = 2^sHkV (15) 

with S ± {k) = J dujS ± (k,uj) the corresponding static structure factor. 

We can establish compressibility sum rules for S' ± (fc,u;) by considering the response of 
the system to an external perturbation H'± = ~ Dk(Pk + P-k) that couples to either the 
even or odd modes. We follow the analysis of Ref. |14| and assume that in the long wave 
length limit the perturbed wave function is locally the same as the unperturbed one with a 
modulated density (even or odd for each mode). We obtain 



lim / du K ; J = 7 ^ r , (16) 



< S ± (k,uj) _L 
fc^o J u 2mv± 

with 

V + = — EI* > V - = — HZ2 > ( 17 ) 



■2 n d 2 €o(n, nJ) 2 n d 2 e (n, n_) 

m On 2 ~ m dn 2 



and eo(n, nJ) the ground state energy per unit volume written as a function of the total 
density n = n a + n& and the density difference n~ = n a — rib. Since in the long wave- 
length limit the sum rules are exhausted by the above quasiparticles [JL5|] [lim fc _> S ,± (/;;, u>) = 
S ± (k)S(cu — co k )}, the sum rules ([H|) and C |X6|) imply that two branches have energies 
uj k = v±k, with the corresponding structure factors given by S ± (k) = k/2mv±. The low 
energy spectrum therefore consists of two linear modes, corresponding to modulations of 
the density with the two species "in-phase" (even mode) and "out-of phase" (odd mode) 
respectively. The odd mode has total density n constant in all space, and in the spinor 
language (where the species index is treated as a spin 1/2) corresponds to a spin density 
wave. We can see qualitatively that the odd modes have lower frequency by perturbing 
around low values of the inter-species interaction. For U x = the even and odd modes are 
degenerate. If we turn on U x the odd mode will have lower expectation value (U x ) since 
locally (p a (x)p&(x))+ > (p a (x)p;,(x))_. This is reasonable as long as the repulsive compo- 
nent of U x is dominant. If we now increase U x we will reach the instability discussed above: 
when U x = U a the ground state is multiply degenerate with all the polarized wave functions 



|*o[(iV - M)/2, (JV + M)/2]) (M = -N, ■ ■ ■ , N) having the same energy. This implies that 
V- —>■ 0, since eo(n, n_) = eo(«.,0). Therefore the instability towards a polarized state is 
signaled by a vanishing velocity of the odd modes, or, equivalently, by a divergence of the 
"compressibility" /t_ oc [d 2 eo{n,n_)/dn 2 _]~ 1 associated with changes in species polarization 
at constant particle number. 

Finally, we note that our proof is not valid in the presence of a field coupling the two 
species of the form lZfd 3 x (ip\(~x)ipb{~x) + h.c). This term corresponds (in the spinor lan- 
guage) to a magnetic field tilted with respect to the direction of the otherwise fully polarized 
system. If the system is prepared in the unstable regime with fixed N a and Nb and no mech- 
anism of interconversion is allowed, our proof indicates that the system will separate into 
two phases. This kind of phase separation was discussed within the Bogoliubov model by 



Nepomnyaschii [Tj|. An infinitesimal coupling with a heat bath that allows the system to 
equilibrate will equate the chemical potentials of both species and induce a transition to a 
state of one species-only The precise kinetics and the time scale to reach equilibrium are 
beyond the scope of our paper. 

When TZ 7^ the solutions will not be completely polarized even in the case where 
jjx > jja^ j.^ ^at the ground state will be "rotated" with respect to the quantization 
axis in species space. The detailed mean field analysis of this problem is the subject of a 



forthcoming paper [17 
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